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P  SPECTRAL  VARIATIONS  OF  THE  GASfiUGGY  EXPLOSION 
n 

AT  INTERMEDIATE  DISTANCE  RANGES 
BY 

W.  H.  K.  Lee  and  R.  D.  Borcherdt 
National  Center  for  Earthquake  Research 
U.S.  Geological  Survey,  Menlo  Park,  California 

INTRODUCTION 

The  Project  GASBUGGY  nuclear  explosive  of  26-kHotons  design 

yield  was  detonated  on  Sunday,  December  10,  1967,  !230t06.!36  MST,  , 

—  /  '  * 

at  latitude  36®40M0.4''  North  and  longitude  Ie7‘"f2‘-30r3''  West,  about 
55  air  miles  east  of  Farmington,  New  Mexico.  The  U.S.  GeologlOal  Survey 
recorded  seismic  waves  generated  by  the  explosion  along  five  lines 
radiating  from  the  shot  site,  us>ng  truck-mounted  seismic-refraction 
recording  systems*(descrlbed-by ■Warrlckr'fit  al r  ( 1961 ),  plus  other 
units.  The  unit  Impulse  response  for  ground  velocity  of  the  record¬ 
ing  systems  Is  approximately  flat  from  2  cps  to  above  15  cps.  A 
preliminary  data  report  on  traveltime  and  amplitude  variations  with 
azimuth  and  distance  was  prepared  by  Warren  and  'Jackson  (1968). 

The  present  paper  Is  submitted  as  a  data-ana lysis  report  whose 
primary  purpose  Is  to  study  variations  of  spectra  as  a  function 
of  distance  and  azimuth.  Seismic  records  were  digitized  from  19  stations 


in  the  distance  range  300  to  500  kilometers  from  the  shotpoint  (Figure''! 
and  Table  I  )-  Time-series  analysis  was  then  appi ied  to  the  portion  of 
seismic  records  containing  essentially  the  phase  (the  first  3  seconds 
after  the  first  arrival).  The  results  Indicated  significant  variations 
of  spectra  from  station  to  station.  Following  a  brief  summary  of 
the  method  of  data  analysis,  the  results  will  be  discussed  in  detail. 


DATA  ANALYSIS 

The  first  3  seconds  of  the  seismic  records  after  the  first 
arrival  were  digitized  at  a  sampling  rate  of  50  per  second,  and 
calibrated  such  that  the  data  represent  absolute  ground  velocity  in 
p/sec  normalized  with  zero  mean  value.  The  usual  method  of  time- 
series  analysis  was  then  applied  as  described  by  Blackman  and  Tukey 
(1958).  For  a  time  series  X(ti),  X(t2),  ...»  X(t^)  with  zero  mean, 
thd  autocovariance  function  at  lag  5  is  computed  by 


R<J)  Xft,) 

for  j  =  0,  1,2,  ...,  m;  andc=  jAt',  where  At  =  tj  - 

The  raw  estimate  of  the  spectral  density  at  frequency  fj^  Is 
Obtained  by 


(I ) 


where 


P(f,) 


2  At 

TT 


rn 

I  e(J)  R(J)  cos(2tt  jAt) 

j  =  0 


0  <  J  <  m 
j.  =  0,m 


(2) 


for  f 


k  =  0,  i. 


m. 
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Figure  I.  Location  map  for  the  GASBUGGY  nuclear  explosion 
and  recording  stations. 


3 


Table  I 


Station  Locations  and  Distances,  Spectral  Energy 
for  the  Phase,  and  Q 


D 1 rect Ion 

Name 

Locat 

Latitude, 

N 

ion 

Longitude, 

W 

Dl stance, 

kn 

Spectral 
Energy, 
(u/sec) 2 

Q 

EAST 

Roy,  N.  M. 

36®03» 

103*56’ 

302 

0.0320 

Rosebud,  N.  M. 

35®52’ 

103*21  ’ 

358 

0.0081 

160 

Romero,  Tex. 

35®43» 

102*56’ 

399 

0.0010 

1  10 

Channing,  Tex. 

35°42» 

102*23’ 

447 

0.0023 

220 

Amari ! lo,  Tex. 

35®3i  ' 

101*49’ 

502 

0.0014 

250 

SOUTH 

Socorro,  N.  M. 

33®58» 

106*57’ 

301 

0.0280 

Mesquite,  N.  M. 

32°I2» 

106*36’ 

499 

0.0012 

250 

SOUTHWEST 

Holbrook,  Ariz. 

35®06’ 

1 10*03’ 

31  1 

0.8400 

Winslow,  ArIz. 

35*04 ' 

110*39’ 

359 

0.1  100 

90 

Chevelon,  Ariz. 

34*37 » 

110*53’ 

404 

0.1  100 

170 

Sunflower,  Ariz. 

33*51  ’ 

1 1 1*26’ 

496 

0.5100 

1200 

WEST 

Shonto,  Ariz. 

36*36’ 

1 10*36’ 

304 

0.3700 

Kaibito,  Ariz. 

36*33’ 

1 10*07’ 

351 

0.1400 

190 

Kanab,  Utah 

37*01 ’ 

112*49’ 

502 

0.0350 

350 

NORTH 

Cl Imax,  Colo. 

39*20’ 

106*13’ 

307 

1  .3000 

Dillon,  Colo. 

39*44’ 

106*03’ 

351 

0.2600 

1  10 

Granby,  Colo. 

40*03’ 

105*56’ 

391 

0.1200 

140 

Timber  Creek,  Colo. 

40*23’ 

105*51 ’ 

428 

0.0013 

70 

Tie  Siding,  Wyo. 

41*02’ 

I05*33‘ 

505 

0.0023 

120 
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The  raw  spectral  density  estimates  are  then  smoothed  by  ’’hamming": 

P(fo)  =  0.54  P(fo)  +  0.46  P(fi) 

P(f^)  =  0.23  P(f^.i)  +  0.54  P(f|^)  +  0.23  0  <  k  <  m 

P(f  )  =  0.54  P(f  )  +  0.46  P(f  ).  (3) 

m  m  m-i 

Spectra  of  ground  velocity  for  the  P^  phase  recorded  along  five 
lines  of  stations  are  shown  In  Figures  2  through  6.  In  each  of  these 
figures,  the  spectral  density  of  ground  velocity  for  the  P^  phase  Is 
plotted  versus  frequency.  Stations  In  the  Same  direction  are  plotted 
In  the  same  figure  to  show  the  spectral  variations  as  a  function  of 
distance.  Dominant  peaks  In  the  spectra  fall  mostly  between  3  to 
5  cpS.  In  general,  the  maximum  spectral  density  decreases  as  the 
distance  Increases  from  the  shotpolnt.  A  notable  exception  Is  station 

o 

Sunflower  at  500  km  where  the  maximum  spectral  density  Is  almos'^  the 
same  as  for  station  Holbrook  at  300  km. 

Spectra  of  ground  velocity  for  the  P^  phase  recorded  at  distances 
of  300  and  500  kilometers  from  the  shotpolnt  are  plotted  In  Figures  7 
and  8  to  show  the  spectral  variations  as  a  function  of  azimuth.  At 
300  km  the  maximum  spectral  density  for  the  north,  west,  and  southwest 
Is  larger  by  a  factor  of  about  30  than  that  for  the  east  and  south. 

At  500  km  the  maximum  spectral  density  Is  highest  toward  the  southwest, 
decreases  by  a  factor  of  about  20  toward  the  west,  and  decreases  by 
an  additional  factor  of  about  10  toward  the  north,  east,  and  south. 

Spectra  of  all  the  stations  show  the  spectral  energy  Is  concen¬ 
trated  between  2  and  8  cps  and  becomes  Insignificant  beyond  10  cps. 
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FREQUENCY  (CPS) 


Figure  2.  spectral  variations  along  the  east  line. 
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Figure  4.  Pp  spectral  variations  along  the  southwest  line. 


SPECTRAL  DENSITY  OF  GROUND  VELOCITY,  (yn/SEC)^/CPS 


FREQUENCY  (CPS) 

Figure  5.  spectral  variations  along  the  west  line. 
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Figure  6.  spectral  variations  along  the  north  line. 
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Figure  7.  Pp  spectral  variations  at  300  km  from  the  shotpoint. 


Figure  8.  spectral  variations  at  500  km  from  the  shotpoint. 
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The  spectral  energy  of  ground  velocity  can  therefore  be  obtained  by 
integrating  the  spectral  density  from  2  to  10  cps: 
rIO 

E  =  P  (f)  df 

h. 

The  results  are  tabulated  in  Table  I.  Spectral  energy  is  plotted 
against  distance  in  Figure  9,  and  a  contour  map  of  iso-spectral 
energy  is  shown  in  Figure  10.  These  figures  show  that  the  spectral 
energy  Is  about  2  orders  of  magnitude  higher  for  the  southwest  and 
west  than  for  the  east  and  south.  The  spectral  energy  along  the 
north  line  is  similar  to  the  southwest  and  west  lines  for  the 
distance  interval  300  to  400  km,  but  Is  similar  to  the  east  and 
south  lines  from  400  to  500  km. 

Attenuation  of  the  phase  can  be  studied  by  calculating  the 
Q  values  between  pairs  of  stations.  The  dimensionless  quantity,  Q, 
IS  the  reciprocal  of  the  specific  attenuation  factor  (Knopoff,  1964) 
and  is  defined  by  Gutenberg  (1958)  as: 

^  kTV  ’  (5) 

where 

k  =  the  coefficient  of  absorption 
T  =  period  of  the  seismic  wave 
V  =  seismic  veiocity. 

The  absorption  coefficient  k  is  defined  by 
^  ■  ^2) 

A2  '  (6) 
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SPECTRAL  ENERGY  OF  GROUND  VELOCITY 
from  2  to  10  CPS  ,{)U/SEC)^ 


Figure  9.  Spectral  energy  of  ground  velocity  from  2  to  10  cps 
for  the  phase  as  a  function  of  distance  along  the 
five  lines  of  recording. 
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where 


A1/A2  =  ratio  of  amplitude  of  seismic  wave  at  location  I  to 
location  2 

Di  =  D2  =  distance  between  location  I  and  location  2. 

For  rough  estimates  of  Q,  we  use  the  spectral  energy  of  ground  velocity 
for  the  phase  from  2  to  1 0  cps  so  that: 


T  =  c -  “0.2  sec  , 

5  cps 


and 


V  8  km/sec 


(7) 


The  results  are  tabulated  in  Table  I.  Values  of  Q  range  from  about 
70  to  350  and  agree  with  the  general  estimate  of  60  to  450  made  by 
Anderson  (1966)  for  the  lower  crust  and  upper  mantle.  The  only 
exception  is  the  0  value  (1200)  between  Holbrook  and  Sunflower.  This 
may  be  due  to  local  conditions  at  Sunflower. 

CONCLUSIONS 

Studies  on  the  P  phase  recorded  at  19  stations  for  the  GASBUGGY 
n 

nuclear  explosion  show  the  following: 

(I)  Spectra  of  ground  velocity  differs  from  station  to  station. 
The  gross  characteristics  are  similar,  but  the  maximum  spectral  density 
may  differ  by  3  orders  of  magnitude.  Dominant  peaks  in  the  spectra 
are  mostly  between  3  and  5  cps. 
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(2)  There  is  a  general  decrease  in  the  spectral  energy  as 
distance  increases  from  the  shotpoint.  Spectral  energy  is  strongly 
dependent  on  the  azimuth,  being  two  orders  of  magnitude  higher  toward 
the  southwest  and  north  than  toward  the  east  and  south. 

(3)  A  contour  map  of  Iso-spectral  energy  indicates  that  atienua- 
tlons  of  seismic  waves  is  higher  in  the  east  and  south  than  in  the 
west  and  southwest.  This  result  is  consistent  with  the  heat-flow 
measurements  showing  higher  heat  flow  toward  the  southeast  than  toward 
the  northwest  (John  Sass,  personal  communication). 

(4)  The  reciprocal  of  specific  attenuation  factor,  Q,  was 
estimated  from  the  spectral  energy  ratios  between  pairs  of  stations. 
Values  of  Q  (70  to  350)  are  consistent  with  the  general  range  estimated 
for  the  lower  crust  and  upper  mantle. 
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